SPH simulations and experimental investigation of water flow through a Venturi meter of rectangular cross-section

The flow of water through a horizontal small-scale Venturi tube of rectangular cross-section is simulated using a modified version of the open-source code DualSPHysics, which is based on Smoothed Particle Hydrodynamics (SPH) methods. Water is simulated using the Murnaghan-Tait equation of state so that weak compressibility is allowed. The hydrodynamics is coupled to a Large-Eddy Simulation (LES) turbulence model. The convergence properties of SPH are improved by adopting a C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{2}$$\end{document}2 Wendland function as the interpolation kernel, increased number of neighboring particles and non-reflective open boundary conditions at the outlet of the Venturi tube. The flow structure and differential pressure as well as the mainstream velocity profiles at different stations are compared with calibrated experimental data. A resolution independence test shows that good convergence to the experimental measurements is achieved using four million particles. At this resolution the simulations predict the experimental centerline velocity profile along the Venturi meter for a volumetric flow rate of ten liters per minutes (lpm) with a root-mean-square error of 4.3%. This error grows to 7.1% when the volumetric flow rate increases to 25 lpm. The predicted differential pressure matches the experimental data with errors varying from 1.4% (for 10 lpm) to 6.8% (for 25 lpm). Cross-sectional velocity profiles within the throat and divergent sections differ from the experimental measurements in less than 5.5%. In general, it is shown that the SPH model can provide an efficient and accurate method for recalibrating flow meters at moderately high Reynolds numbers instead of using costly experimental tests.

to their emission into the environment [5][6][7][8][9] .In addition to being one of the most popular and highly efficient gas cleaning devices, Venturi injectors have also found important applications in sprinkler and drip irrigation systems for intensive farming and industrial agriculture 10 .They have also been used in applications of liquid flow metering 11 , mixing of liquid/liquid and gas/liquid systems as, for example, in the injection of ozone in water and in water aeration systems [12][13][14] , and in the enhancement of the chemical reaction of hydrodynamic cavitating flows 15 .Efforts have also focused on studying the instability of cavitating flows in a Venturi meter, whose control can help to limit structural damage 16 .In general, Venturi injectors operate over a wide range of pressures, are simple in design, have no moving parts and are almost free of maintenance.Due to all these advantages and their widespread applications in the industry, a wide range of models of varying complexity have been presented in the literature, which have been frequently validated experimentally.
In addition to the large number of existing experiments, a great deal of effort has also been devoted to the numerical modeling of Venturi flows in order to reproduce the experimental data for laboratory-scale and industrial-size devices.Theoretical studies on the hydrodynamics of Venturi meters range from simple correlations to more complex numerical simulations.For instance, a boundary layer growth model was developed by Gamisans et al. 17 to investigate the pressure drop characteristics of Venturi scrubbers, while the hydrodynamics of a liquid/liquid phase distribution through a Venturi meter in a vertical pipe was examined by Jana et al. 18 .Numerical simulations of flow behavior within a Venturi meter have been performed to determine air injection rates 19 , pressure differences on a bluff body object for water flow 20 and on turbine blades placed within the throat section for airflow 21 , as well as to determine the effects of varying the geometry of the convergent section on pressure drop 22 and study the behavior of two-phase flow and mixing [23][24][25][26][27] .More recently, cavitating flows in a small-sized Venturi tube have been studied by Razali et al. 28 Further numerical simulations have been carried out to measure discharge coefficients and pressure drops in industrial-sized Venturi meters for oxygenfree nitrogen 29 and water flow 30 .The effects of varying the structural parameters of a Venturi tube, including the contraction ratio, the ratio of the throat section length to diameter, the diffusion angle and the inlet and outlet pressure difference, on fluid flow were analyzed by Zhang 31 using the Fluent v6.3 software.Finite element methods have been used for design optimization analysis of a Venturi tube for medium conveying in strengthen grinding processes 32 , while investigation of large over-reading of the Venturi flow rate in the feed water flow control system in the second circuit of a nuclear power plant was recently performed by Wang et al. 9 with the aid of the software package ANSYS Fluent 17.0.
Most early and recent simulations of the Venturi flow have been carried out using traditional mesh-dependent methods and no calculations have been performed to date using the method of Smoothed Particle Hydrodynamics (SPH).In this study, we perform three-dimensional calculations of water flow at Reynolds numbers between 13272 and 33180 through a small-scale Venturi tube of rectangular cross-sectional shape, using a modified version of the open-source code DualSPHysics 33 , which is based on SPH theory for solving the equations of fluid dynamics.Water is simulated using the Murnaghan-Tait equation of state so that a weakly compressible SPH scheme is proposed 34 .The Large-Eddy Simulation (LES) turbulence model is implemented to filter the equations of fluid dynamics.The convergence properties of SPH are improved in three different ways.First, a C 2 Wendland function 35,36 is used as the interpolation kernel.In contrast to most conventional kernels, Wendland functions support large numbers of neighboring particles without becoming unstable and do not allow particle motions on a sub-resolution scale (i.e., within the kernel support).Therefore, they maintain quasi-ordered particle distributions even in highly dynamical test problems.Second, approximate volume partitioning is achieved using a large number of neighbors 37,38 .Third, non-reflecting outlet boundary conditions are implemented by allowing the particles that leave the Venturi diffuser and enter the outflow zone to move according to an outgoing wave equation for the velocity field.This effectively reduces feedback noises from the outlet boundary, improving the accuracy of the solution.The Venturi tube was designed for experimental purposes and the numerically obtained flow behavior is compared with calibrated experimental flow measurements.The motivation for this study is twofold.First, understand the flow measurement device and second, explore the efficiency and accuracy of the proposed SPH method for recalibrating flow meters at moderately high Reynolds numbers.In particular, the flow structure along the convergent, throat and divergent sections of the Venturi as well as the differential pressure and mainstream velocity profiles at different stations along the tube are compared with the experimental data.The results show that the experimental measurements are reproduced with errors less than about 7.1% in the velocity profiles (when measured in units of m s −1 ) and less than ∼ 6.8 % in differential pressure (in units of psi) along the length of the Venturi for the highest flow rates considered.

Experimental flow loop
The experimental flow measurements were performed using a closed flow loop as shown schematically in the upper part of Fig. 1.The Venturi tube is installed horizontally with the differential pressure and has a rectangular cross-sectional area as depicted in the lower part of Fig. 1.It is composed of a convergent section of length 35.45 mm and convergent angle of 30 • , a throat of length 6 mm and a diffuser of length 108.57mm and divergent angle of 10 • .The cross-sectional area at the inlet of the convergent section and outlet of the diffuser is a × b = 6.35 × 25 mm 2 , while the throat section has a cross-sectional area of a × b = 6 × 6.35 mm 2 .The inlet and outlet planes of the Venturi meter have hydraulic diameters, D H = 2ab/(a + b) ≈ 10.13 mm, while the throat has a hydraulic diameter of ≈ 6.17 mm.For this meter, the contraction ratio, calculated as the ratio between the height of the throat and that of the inlet plane, is γ = 6/25 = 0.24.
The Venturi meter is made of transparent acrylic to facilitate flow visualization, while the pipeline upstream and downstream of the Venturi is made of PVC and has a circular cross-section of diameter 25 mm.The procedure undertaken in the experiments is as follows.First, a peripheral pump of 1 HP capacity is switched on to drive water from a 50-liter storage tank open to the atmosphere through an orifice.The flow is controlled by a check valve, which is placed upstream of a flow meter which in turn adjusts the volumetric flow rate at the inlet of the Venturi to be in the range between 10 and 35 liters per minute (lpm).When the circuit is completely filled with water, the pump keeps the circulation going by maintaining a head difference and the desired flow.The dynamic pressure is measured by means of three pressure transducers placed at the inlet (P 1 ), in the middle of the throat (P 2 ), and at the outlet plane (P 3 ) (see Fig. 2).The characteristics of the flow were determined using the particle image velocimetry (PIV) technique, which enables an accurate spatial and temporal resolution of the velocity field.A centerline pressure profile along the length of the Venturi tube is determined by applying Bernoulli's equation for j = 1, 2, • • • , M − 1 , where for j = 1 the value of p j is given by the experimentally measured pressure at point P 1 , the values of v j correspond to the measured mainstream flow velocities along the Venturi meter, and M stands for the number of points on the centerline where measured values of the velocity are known.This procedure is repeated by iterating backward from the value of pressure at point P 3 to check the accuracy in reproducing  the experimentally obtained pressure at point P 2 .The final profile is calculated by averaging the forwardly and backwardly iterated pressures at each available centreline position along the Venturi meter.Images were recorded with a high-speed Phantom Speedsense 9040 camera placed in front of the Venturi tube with an image resolution of 1260 × 1024 pixels and working at 1000 frames per second.The velocity profiles were obtained by seeding the water with polyamide spherical tracer particles of 50 nm diameter each.The Dantec Dynamic Studio software was employed to process the acquired images.The experiments were performed at an ambient temperature of 30 • C for volumetric flow rates between 10 and 25 lpm, corresponding to bulk velocities at the inlet of the Venturi meter v b = 1.05 and 2.625 m s −1 , respec- tively.At such flow rates, the flow through the Venturi meter corresponds to Reynolds numbers, Re= ρv b D H /µ , between 13272 (for Q = 10 lpm) and 33180 (for Q = 25 lpm), where µ = 0.000798 kg m −1 s −1 is the dynamic viscosity of water at 30 • C.

Governing equations
The differential equations describing the flow of water through the Venturi meter are given by the Navier-Stokes equations where ρ is the mass density, v is the velocity vector, p is the pressure, ν is the coefficient of kinematic viscosity, g is the gravitational acceleration and d/dt = ∂/∂t + v • ∇ is the material time derivative.The pressure is defined in terms of the density using the Murnaghan-Tait equation of state 34 where γ = 7 is appropriate for water, p 0 = c 2 0 /γ , ρ 0 = 995.71kg m −3 is the water density at 30 • C and c 0 is the sound speed at the reference density ρ 0 .The kinematic viscosity of water at 30 Here we adopt a weakly compressible approach and the term p 0 governs the relative density fluctuations |ρ − ρ 0 |/ρ 0 ∼ M 2 , where the Mach number M is typically set to 0.01 in order to enforce density fluctuations less than about 1%.This is accomplished by defining the reference sound speed c 0 to be at least ten times higher than the maximum fluid velocity across the flow meter.

Numerical methods
In the SPH framework a LES filtering of Eqs.(1) and ( 2) coupled with a sub-particle scaling technique is implemented to model coherent turbulent structures in the fluid 39 .In this approach, the velocity field is separated into its mean component, ṽ , and its fluctuating part, v ′ , such that v = ṽ + v ′ , where the mean velocity is obtained by means of the density-weighted Favre-filtering where T is a sufficiently large time interval and ρ is the Reynolds-averaged density.After application of the Favre-filtering, Eqs.(1) and (2) become Here, T is the sub-particle stress tensor where is the Favre-filtered strain rate tensor, C I = 0.00066 , ν t = (0.12∇) 2 |S| is the Smagorinsky eddy viscosity, |S| = (2S ij S ij ) 1/2 is the local strain rate, δ ij is the Kronecker delta and is a measure of the finite particle size.
Equations ( 5) and ( 6) are solved numerically using a modified version of DualSPHysics 33 , which is based on SPH methods 40,41 .SPH is a fully Lagrangian, mesh-free scheme for the simulation of complex fluid-flows, where the fluid is represented by a discrete number of moving particles that carry all field information of the system (1) www.nature.com/scientificreports/as, for example, the density and velocity.In DualSPHysics the density of particle a is calculated using Eq. ( 5), which in discrete form is written as where ρ a is the particle-scale density associated to particle a, m b is the mass of neighboring particle b, W ab = W(|x a − x b |, h) is the kernel function, h is the smoothing length and n is the number of neighbors of particle a within the kernel support, i.e., within distances |x a − x b | ≤ h from particle a.The value of n is chosen to be n = 0.12 √ N , where N is the total number of SPH particles filling the computational domain.This prescrip- tion complies with the joint limit N → ∞ , n → ∞ and h → 0 for complete convergence to the continuum 37,38 .The SPH form of Eq. ( 6) is given by where x ab = x a − x b and ε 2 = 0.01h 2 .In the above SPH representations the bar and tilde operators over the mean density and ensemble average velocity vector have been dropped for simplicity.
The particles are moved along the Venturi meter by integrating the equation where the second term on the right-hand side prevents the growth of errors due to anisotropies in the distribution of particle positions 42 .Here β = 0.04 , v max is the maximum velocity, M is the total fluid mass within the computational domain and x 0 , given by is the mean distance between particle a and all other particles, where the summation is taken over all N particles filling the computational domain.The convergence properties of SPH are improved by adopting a C 2 Wendland function as the interpolation kernel 35,36 for q ≤ 1 and zero otherwise, where q = |x − x ′ |/h .The time integration of Eqs. ( 9)-( 11) is performed using the Verlet algorithm provided by DualSPHysics, which is second-order accurate and maintains adequate numerical coupling among Eqs.( 9)-( 11) during the evolution.The use of Eq. ( 3) enforces a weakly compressible SPH scheme.At small scales most fluids are incompressible and an alternative SPH approach that satisfies this condition is the so-called incompressible SPH 43 .This approach has been successfully applied in the simulation of free-surface flows, including Newtonian and non-Newtonian flows, dam break flows and wave overtopping on different coastal structures [44][45][46][47] , among many other applications.The numerical integration of Eqs. ( 9)-( 11) is performed by applying no-slip boundary conditions (i.e., v = 0 ) at the solid walls of the Venturi tube using the method of dynamic boundary particles proposed by Crespo et al. 48,49 .In these simulations, a layer of particles at rest defines the walls of the Venturi and two layers of ghost particles are placed outside the computational domain around the volume occupied by the Venturi meter.Although these particles are updated using Eq. ( 10), they are not allowed to move so they preserve their initial positions.The density of the outer ghost particles is determined from the neighboring fluid particles lying within their support domains.Particle penetration across the solid walls is avoided by allowing the wall particles to exert repulsive forces on nearby fluid particles, which are derived from the source of Eq. (10).
Open boundary conditions are applied at the inlet, L 1 , and outlet, L 5 , planes (see Fig. 2).An inflow zone is placed in front of the inlet plane, where inflow particles are arranged regularly and allowed to flow in as needed with prescribed density and velocity values 50 .Non-reflecting outflow boundary conditions are implemented by defining an outflow zone in front of the outlet (L 5 ) plane such that particles there are moved using an outgoing wave equation of the form where v = (v x , v y , v z ) and the mainstream velocity across the L 5 -plane is assumed to be along the x-axis.A numerically stable SPH representation of this equation can be written as ( 9) where the subscript o is used to denote particles in the outflow zone.Here Since the inlet and outlet mass rates may differ from each other, a particle leaving the outflow zone is temporarily stored in a reservoir buffer.As an inflow particle crosses the L 1 -plane, a particle is removed from the buffer and placed in the upstream side of the inflow zone with a prescribed density and velocity.A more stable scheme with the use of Eq. ( 15) can be obtained by smoothing the mainstream velocity component, v x,o , according to the prescription where a particle o can have neighbors pertaining to the fluid and outflow zones depending on how close it is from the outlet plane.Finally, the position and velocity of outflow particles is obtained by integrating the equation simultaneously with Eq. ( 15) using the same Verlet integrator employed to evolve the hydrodynamics.

Convergence testing and numerical validation
A resolution independence test was performed for the SPH simulations to determine the quality of the numerical solutions and convergence to the experimental data.The experimentally obtained mainstream velocity variation along the Venturi meter for a volumetric flow rate of 10 lpm was chosen as a convergence test for numerical validation purposes.The SPH particles are initially at rest and uniformly distributed in all three coordinate directions so that = x = y = z and the smoothing length is set equal to C √ 3 , with C = 1.15 , for all runs.Figure 3 shows SPH results as the number of particles is increased from N = 254288 to 4054963 particles.The numerical profiles show a tendency to globally converge to the experimental data (open circles) as the spatial resolution is improved.The last two columns of Table 1 list the maximum velocity at the exit of the throat section ( x = 0.115 m) and the root-mean-square errors (RMSEs) between the numerical and experimental profiles at different resolutions, respectively.In terms of this metric, the error between the SPH and the experimental data decays from ≈ 0.41 m s −1 for N = 254288 particles (with more significant deviations in the Venturi diffuser) to ≈ 0.04 m s −1 for N = 4054963 particles (see Table 1).At the highest resolution the numerical profile fits the experimental one with only an error of ≈ 4.3 %.The position and magnitude of the experimental peak velocity ( ≈ 2.887 m s −1 ) at x = 0.115 m is very well reproduced for N ≥ 2015768 particles, with values of ≈ 2.89 m s −1 (for N = 2015768 ) and ≈ 2.92 m s −1 (for N = 4054963 ).The higher dispersion of the SPH velocity profiles along the diffuser ( x < 0.075 cm) for ≤ 1 million particles occurs because near the pipe walls the SPH velocity values are larger than the real ones, resulting in more pronounced parabolic profiles than it should be.As the number of particles is increased, the flow velocity near the walls becomes closer to the experimental measurements and the dispersion is drastically reduced.( 15)  2 list the RMSEs between the experimental and SPH velocity and pressure profiles, respectively.These errors increase with increasing volumetric flow rate.For the velocity, the SPH simulations reproduce the experimental data with ≈ 4.3 % for Q = 10 lpm and ≈ 7.1 % for Q = 25 lpm.When measured in units of psi, the pressure profiles are reproduced with RMSEs that vary between ≈ 1.4 % for Q = 10 lpm to about 6.8% for Q = 25 lpm.The RMSE between the experimentally measured and the numerically obtained differential pressure also increases with the volumetric flow rate.
The experimentally observed position of the maximum velocity and pressure at the exit of the throat section ( x ≈ 0.115 m) is very well reproduced by the SPH simulations.The experimental and SPH maximum veloci- ties in units of m s −1 are ( v max,exp ,v max,SPH )=(2.89,2.92)for Q = 10 lpm, (6.12,6.14)for Q = 15 lpm, (7.86,7.75)for Q = 20 lpm, and (9.30,9.53)for Q = 25 lpm.These figures show that the strength of the peak velocities is fairly well reproduced by the SPH simulations.The largest differences are observed for Q = 20 and 25 lpm.The numerically obtained maximum velocities as compared to the experimental data are depicted in Fig. 5 for all volumetric flow rates considered.In units of psi, the experimental and SPH minimum pressure values at the exit of the throat section are: ( p min,exp ,p min,SPH )=(−0.70,−0.72 ) for Q = 10 lpm, ( −1.63,−1.57) for Q = 15 lpm, ( −2.32,−2.16 ) for Q = 20 lpm, and ( −2.75,−2.88 ) for Q = 25 lpm.
The mainstream velocity profiles in the horizontal (z) plane (see Fig. 1), with z ⋆ = z/(6.35mm), as compared with the experimental measurements are shown in Fig. 6     were repeated several times for reproducibility and so they are considered to be accurate enough to provide a useful benchmark test case for checking the accuracy of the SPH flow prediction method.In all cases the SPH profiles closely match the experimental data with RMSE deviations ranging from ≈ 2.1 % at L 3 for Q = 10 lpm to about 5.5% at L 4 for Q = 20 lpm (see Table 3).As the core of the vena contracta expands through the diffuser, the velocity towards the pipe walls decays less steeply consistently with the flow there approaching a parabolic shape.Table 3 lists the RMSE deviations between the experimental and SPH cross-sectional mainstream velocity profiles for all volumetric flow rates considered.The errors increase with increasing values of Q.At pipe station L 3 , within the throat section, the profiles are almost flat with very sharp decays to zero close to the pipe walls.As Q is increased the extension of the vena contracta also increases and for Q = 25 lpm it becomes as long as the  divergent section.This explains the almost flat profile in the central core with less steep decays to zero towards the pipe walls for Q ≥ 15 lpm.In particular, for Q = 10 lpm the vena contracta downstream of the throat exit shortens and the flow becomes almost flat, though at a lower velocity compared to station L 3 .Close to the exit of the diffuser (station L 5 ), the vena contracta weakens and the flow try to revert to laminar.However, a smaller fraction of the central core still remains approximately flat.Therefore, as the vena contracta completely diffuses outside the Venturi tube the profile gradually approaches a parabolic shape.Figure 7 shows a direct comparison between the experimental and the SPH calculated mainstream velocity field along the Venturi meter for (a) 15 lpm and (b) 25 lpm.The velocity fields are qualitatively similar with small differences close to the inlet and outlet of the Venturi tube.

Conclusions
In this study a modified version of the open-source code DualSPHysics, which relies on the method of Smoothed Particle Hydrodynamics (SPH), has been used to simulate the flow of water at moderately high Reynolds numbers through a horizontal, small-scale Venturi meter of rectangular cross-sectional shape.The geometry and dimensions of the Venturi tube as well as the initial flow parameters used for the numerical simulations correspond to the experimentally defined values.
As the flow rate is increased within the allowed experimental range, the SPH simulations predict mainstream velocity profiles and differential pressures along the Venturi tube that are in good qualitative agreement with the experimental data.The deviations between the SPH and the experimental mainstream velocity profiles vary from 4.3% for a volumetric flow rate (Q) of 10 liters per minute (lpm) to 7.1% for Q = 25 lpm, while the SPH differential pressures differ from the experimental data by about 1.4% (for Q = 10 lpm) and 6.8% (for Q = 25 lpm).The cross-sectional velocity profiles within the throat and divergent sections are always less than 5.5%.In general, the largest deviations from the experimental data occur for Q = 25 lpm.At higher flow rates, additional experiments have shown that cavitating flows may form within the throat section due to much stronger pressure losses.As a further step towards improving the quality of the SPH simulations for calibrating laboratory-scale and industrial-size flow metres, the results from this study will be extended to simulate higher flow rates and include thermal effects through the use of an internal energy equation at much higher spatial resolutions to allow for direct numerical simulations.The present results have applications in deep water culture (DWC) systems 51 .Section varying devices with square or rectangular cross-section are common implementations in controlled irrigation systems within the schemes of precision farming as in open-loop techniques and DWC, due to the square geometry of flooded receptacles for plants.As a low-impact tool, the Venturi meter studied here can be

1 .
Figure 1.Schematic of the flow test loop used for the experiments (top) and three-dimensional view and description of the small-sized Venturi flow meter (bottom).

Figure 2 .
Figure2.Schematic side view of the Venturi meter.P 1 , P 2 and P 3 mark the points where pressure values were measured, while L 1 , L 2 , L 3 , L 4 and L 5 mark the stations along the Venturi where mainstream cross-sectional velocity profiles were measured experimentally.

o dt = v o , Figure 3 .
Figure 3. Centerline velocity profiles along the Venturi meter for Q = 10 lpm at different spatial resolutions as compared to the experimental data.Asymptotic convergence to the experimental data is already observed for N ≥ 2015768 SPH particles.For reference x = 0.15 m and x = 0 m mark the positions of the entrance (L 1 ) and exit (L 5 ) planes of the Venturi meter, respectively.
Figure 4a,b show the centerline velocity and pressure profiles along the Venturi meter for Q = 10 , 15, 20 and 25 lpm.In each figure the position x = 0.25 m marks the entrance plane (L 1 ), while x = 0 m marks the exit plane (L 5 ) of the Venturi tube.The last two columns of Table2list the RMSEs between the experimental and SPH velocity and pressure profiles, respectively.These errors increase with increasing volumetric flow rate.For the velocity, the SPH simulations reproduce the experimental data with ≈ 4.3 % for Q = 10 lpm and ≈ 7.1 % for Q = 25 lpm.When measured in units of psi, the pressure profiles are reproduced with RMSEs that vary between ≈ 1.4 % for Q = 10 lpm to about 6.8% for Q = 25 lpm.The RMSE between the experimentally measured and the numerically obtained differential pressure also increases with the volumetric flow rate.The experimentally observed position of the maximum velocity and pressure at the exit of the throat section ( x ≈ 0.115 m) is very well reproduced by the SPH simulations.The experimental and SPH maximum veloci- ties in units of m s −1 are ( v max,exp ,v max,SPH )=(2.89,2.92)for Q = 10 lpm, (6.12,6.14)for Q = 15 lpm, (7.86,7.75)for Q = 20 lpm, and (9.30,9.53)for Q = 25 lpm.These figures show that the strength of the peak velocities is fairly well reproduced by the SPH simulations.The largest differences are observed for Q = 20 and 25 lpm.The numerically obtained maximum velocities as compared to the experimental data are depicted in Fig.5for all volumetric flow rates considered.In units of psi, the experimental and SPH minimum pressure values at the exit of the throat section are: ( p min,exp ,p min,SPH )=(−0.70,−0.72 ) for Q = 10 lpm, ( −1.63,−1.57) for Q = 15 lpm, ( −2.32,−2.16 ) for Q = 20 lpm, and ( −2.75,−2.88 ) for Q = 25 lpm.The mainstream velocity profiles in the horizontal (z) plane (see Fig.1), with z ⋆ = z/(6.35mm), as compared with the experimental measurements are shown in Fig.6at pipe stations (a) L 3 , (b) L 4 and (c) L 5 (see Fig. 2), respectively.Each figure depicts the profiles for Q = 10 , 15, 20 and 25 lpm.The experimental measurements

Figure 4 .
Figure 4. (a) Centerline mainstream velocity and (b) pressure profiles along the Venturi meter for varying volumetric flow rates between Q = 10 and 25 lpm as compared with the experimental data.

Figure 5 .
Figure 5. Experimentally measured maximum flow velocities at the exit of the throat section (solid line) as compared with the numerically obtained values (blank squares).

Figure 6 .
Figure 6.Numerical horizontal mainstream velocity profiles for varying volumetric flow rates between Q = 10 and 25 lpm as compared with the experimental data at pipe stations: (a) L 3 within the throat section, (b) L 4 in the middle of the divergent section and (c) L 5 at the exit plane of the divergent section.

Figure 7 .
Figure 7.Comparison between the experimental and SPH mainstream velocity maps along the Venturi meter for (a) 15 lpm and (b) 25 lpm.

Table 1 .
Spatial resolution parameters, maximum velocity and RMSEs for the convergence test.

Table 2 .
Flow characteristics and errors of centerline mainstream velocity and pressure profiles along the Venturi meter.

Table 3 .
RMSE deviations between the experimental and SPH cross-sectional mainstream velocity profiles.